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PREFACE 

The  work  reported  herein  was  conducted  by  the  Arnold  Engineering  Development 
Center  (AEDC),  Air  Force  Systems  Command  (AFSC).  The  results  of  the  research  were 
obtained  by  ARO,  Inc.,  AEDC  Group  (a  Sverdrup  Corporation  Company),  operating 
contractor  for  AEDC,  AFSC,  Arnold  Air  Force  Station,  Tennessee.  The  AEDC  project 
manager  was  Mr.  Alexander  F.  Money.  The  research  associated  with  the  analysis  of  test  data 
and  with  the  development  of  a  motion  simulation  computer  program  was  conducted  under 
ARO  Project  Number  P32F-30C.  Testing  was  conducted  from  March  27  through  April  5, 
1979,  under  ARO  Project  Number  P41C-A3.  The  manuscript  was  submitted  for  publication 
on  January  25,  1980. 
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1.0  INTRODUCTION 

The  flight  envelope  of  current  fighter  aircraft  includes  maneuvers  at  high  angles  of  attack 
at  transonic  speeds.  The  prediction  of  flight  characteristics  for  these  conditions  must  rely  on 
flight  test  or  wind  Tunnel  data  that  may  include  nonlinearities  and,  in  some  cases,  hysteresis, 
which  is  defined  herein  as  the  condition  of  dependence  of  data  values  on  the  past  history  of 
the  model  attitude.  The  presence  of  hysteresis  was  once  attributed  to  model  asymmetries  or 
discrepancies  that  were  not  present  on  the  full-scale  aircraft;  however,  this  has  since  been 
disproved  (see  Ref.  1). 

Hysteresis  in  static  aerodynamic  data  obtained  in  wind  tunnel  tests  of  several  fighter 
designs  has  been  documented  in  Ref.  2.  The  longitudinal  and  lateral/directional  data 
exhibited  hysteresis  loops  as  well  as  large  nonzero  values  of  coefficients  at  zero  sideslip 
angles.  These  effects  were  noted  in  the  mid-to-high  angle-of-attack  range  (from  15  to  30 
deg),  primarily  in  the  transonic  Mach  number  range.  Two  of  the  models  in  Ref.  2  —  a 
0.05-scale  model  of  the  F-4,  and  a  0.0417-scale  F-l  1 1  model  —  were  chosen  for  further  tests 
to  study  the  hysteresis  phenomena;  these  tests  were  conducted  in  the  Aerodynamic  Wind 
Tunnel  (4T)  of  the  Arnold  Engineering  Development  Center  (AEDC)  Propulsion  Wind 
Tunnel  (PWT)  facility  and  are  documented  in  Ref.  3.  To  investigate  the  cause  of  hysteresis, 
static  longitudinal  and  lateral/directional  force  and  moment  data  and  wing  pressure  data 
were  obtained  for  the  F-4  model  (without  pylons  or  external  stores)  and  with  various 
simulated  wing  leading-edge  slats.  These  data  were  obtained  at  Mach  numbers  from  0.7  to 
0.95.  Angle  of  attack  was  varied  from  -2  to  24  deg  at  zero  sideslip  angle,  and  sideslip  angle 
was  varied  from  -12  to  12  deg  at  angles  of  attack  of  5,  10,  15,  and  20  deg.  Selected  data  from 
the  F-4  model  test  are  presented  and  discussed. 

In  addition,  computer-generated  six-degree-of-freedom  motion  simulation  studies  were 
used  to  assess  the  effects  of  rolling-moment  hysteresis  on  the  prediction  of  aircraft  motion. 
Simulations  of  various  flight  maneuvers  were  conducted  both  with  and  without  hysteresis  in 
the  aerodynamic  data.  Several  “logics”  were  developed  for  applying  the  wind  tunnel  data  in 
the  simulation  program. 

2.0  WIND  TUNNEL  TEST  PROGRAM  APPARATUS 
2.1  TEST  FACILITY  AND  MODEL  SUPPORT  SYSTEM 

The  AEDC  Aerodynamic  Wind  Tunnel  (4T)  is  a  closed-loop,  continuous-flow,  variable- 
density  tunnel  in  which  Mach  number  can  be  varied  from  0.1  to  1.3  and  set  at  discrete  Mach 
numbers  of  1.6  and  1.96  by  placing  nozzle  inserts  over  the  permanent  sonic  nozzle.  At  all 
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Mach  numbers,  the  stagnation  pressure  can  be  varied  from  300  to  3,700  psfa.  The  test 
section  measures  4  ft  square  by  12.5  ft  long,  and  it  has  perforated  walls  whose  porosity  can 
be  varied  from  0.5  to  10  percent  open.  It  is  completely  enclosed  in  a  plenum  chamber  from 
which  part  of  the  tunnel  airflow  can  be  evacuated  through  the  perforated  walls.  The  model 
support  system  consists  of  a  sector  and  sting  attachment  which  has  a  pitch  angle  capability 
of  from  -8  to  27  deg  with  respect  to  the  tunnel  centerline  and  a  roll  capability  of  from  -180  to 
180  deg  about  the  sting  centerline.  (A  more  complete  description  of  the  tunnel  may  be  found 
in  Ref.  4.) 

2.2  TEST  ARTICLE 

The  test  article  was  a  0.05-scale  model  of  the  F-4C  aircraft.  Details  of  the  F-4C  model  are 
shown  in  Fig.  la.  The  left  wing  of  the  model  was  constructed  with  static  pressure  orifices 
located  at  two  spanwise  stations  on  the  upper  surface  as  shown  in  Fig.  lb.  The  simulated 
leading-edge  slats  that  were  tested  with  the  basic  wing  are  shown  in  Fig.  lc.  The  outboard 
slat,  Si,  extended  from  the  wing  leading-edge  notch  {BL  8.00)  to  BL  11.03.  A  mid  slat,  S2, 
extended  from  BL  4.75  to  BL  8.00.  This  slat  was  tested  in  combination  with  S|.  A  slat 
designated  as  S3  extended  from  the  wing/fuselage  junction  to  the  wing  leading-edge  notch. 
The  leading-edge  slats  were  attached  to  the  basic  wing;  no  attempt  was  made  to  simulate  the 
actual  wing  leading  edge  with  slats  extended.  The  model/slat  configurations  are  identified  in 
Table  1.  The  model  stabilator  was  held  constant  at  0  deg  with  respect  to  a  waterline;  the 
model  had  flow-through  inlet  ducts.  The  model  installation  in  Tunnel  (4T)  is  shown  in  Fig. 
2.  Free  boundary-layer  transition  was  allowed  on  all  model  components  throughout  the  test. 

2.3  INSTRUMENTATION 

A  six-component,  internal  strain-gage  balance  was  used  to  measure  the  aerodynamic 
forces  and  moments  on  the  F-4C  model.  Model  pressure  measurements  were  made  with  use 
of  differential  pressure  transducers  of  the  Tunnel  (4T)  pressure  system. 

3.0  TEST  DESCRIPTION 
3.1  TEST  CONDITIONS  AND  PROCEDURES 

Static  aerodynamic  data  were  obtained  for  all  configurations  at  Mach  numbers  from  0.7 
to  0.95  at  a  constant  total  pressure  of  1,200  psfa.  The  nominal  test  conditions  were 
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M„ 


Re  x 


,  per  foot 


0.7 

2.18 

0.8 

■  2.32 

0.85 

2.37 

0.9 

2.42 

0.95 

2.51 

The  data  were  obtained  by  varying  model  pitch  and  roll  angle  and  by  utilizing  computer 
facilities  to  set  the  model  pitch  and  roll  angles  to  give  the  required  values  of  angle  of  attack 
and  sideslip  angle  at  a  fixed  test  condition.  Extreme  care  was  taken  throughout  the  test  to 
monitor  and  document  the  model  attitude  history  in  order  to  correctly  interpret  the  results 
when  hysteresis  was  present  in  the  data. 

Pitch  hysteresis  data  are  static  aerodynamic  data  that  are  obtained  while  varying  the 
angle  of  attack  at  a  constant  angle  of  sideslip,  and  yaw  hysteresis  data  are  static  aerodynamic 
data  that  are  obtained  while  varying  the  angle  of  sideslip  at  a  constant  angle  of  attack.  For 
pitch  polars,  data  were  obtained  as  the  model  angle  of  attack  was  varied  from  the  minimum 
to  the  maximum  angle  of  attack  and  then  returned  to  the  minimum  angle.  Yaw  polars,  with 
the  exception  of  model  attitude  history  studies,  were  obtained  by  (1)  starting  the  model  at 
zero  deg  angle  of  attack  and  sideslip  angle,  (2)  pitching  the  model  to  the  desired  angle  of 
attack  at  zero  sideslip  angle,  (3)  moving  the  model  to  the  maximum  negative  sideslip  angle, 
(4)  moving  to  the  maximum  positive  sideslip  angle,  and  (5)  returning  to  the  maximum 
negative  sideslip  angle.  The  arrows  on  some  of  the  figures  indicate  direction  of  model 
movement. 

3.2  DATA  REDUCTION  AND  CORRECTIONS 

All  force  and  moment  data  presented  in  this  report  are  measured  body  axis  coefficients. 
Moments  were  referenced  to  FS  16.217,  WL  1.55,  and  BL  0  (see  Fig.  la). 

The  angle  of  attack  and  angle  of  sideslip  were  corrected  for  sting  and  balance  deflections 
caused  by  aerodynamic  loads.  Corrections  for  components  of  weight,  normally  termed  static 
tares,  were  also  applied  to  the  data.  The  model  was  tested  both  upright  and  inverted  to 
provide  the  data  to  correct  for  the  tunnel  flow  angularity  in  the  pitch  plane. 
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4.0  TEST  RESULTS  AND  DISCUSSION 

4.1  HYSTERESIS  CHARACTERISTICS  IN  PITCH 

Figure  3  presents  the  longitudinal  characteristics  of  the  clean  F-4C  model,  showing  the 
presence  of  aerodynamic  hysteresis.  Hysteresis  is  present  for  normal-force,  pitching- 
moment,  axial-force,  and  rolling-moment  coefficients  over  a  wide  range  of  angles  of  attack 
and  Mach  numbers  (see  Figs.  3a,  b,  c,  and  f).  The  variations  of  wing  upper  surface  pressure 
distribution  with  angle  of  attack  (with  increasing  a)  are  presented  in  Fig.  4  at  =  0.7. 
Figure  5  presents  the  effect  of  model  pitch  direction  on  wing  pressure  distributions  (clean 
F-4C)  for  a  =  15  deg  (with  the  exception  of  =  0.7,  where  a  =  14  deg).  Wing  pressure 
distributions  of  the  outboard  section  (BL  9.76)  were  dependent  on  the  direction  of  model 
movement  for  all  subsonic  Mach  numbers  (M^  -  0.7  through  0.95).  The  flow  separation 
and  reattachment  characteristics  were  also  dependent  on  Mach  number.  At  the  higher 
subsonic  Mach  numbers  (Ma  =  0.85  to  0.95),  the  flow  on  the  outboard  section  of  the  wing 
upper  surface  is  completely  separated  at  a  =  15  deg  with  decreasing  angle  of  attack  as 
shown  in  Figs.  5c  through  e.  The  flow  separation  and  reattachment  were  not  symmetrical,  as 
indicated  by  the  rolling-moment  coefficient  (see  Fig.  3f). 

4.2  HYSTERESIS  CHARACTERISTICS  IN  YAW 

The  effect  of  yaw  angle  on  the  lateral/directional  characteristics  of  the  clean  F-4C  at 
M.  =  0.9  is  presented  in  Fig.  6  for  various  angles  of  attack.  Hysteresis  was  predominant  in 
the  rolling-moment  coefficient  at  a  15-deg  angle  of  attack,  and  the  value  of  the  rolling- 
moment  coefficient  was  dependent  on  the  previous  attitude  history  of  the  model  (see  Fig. 
6c). 

The  effect  of  Mach  number  on  the  hysteresis  characteristics  of  the  rolling-moment 
coefficient  is  shown  in  Fig.  7  for  a  model  angle  of  attack  of  15  deg.  Hysteresis  in  the  rolling- 
moment  coefficient  occurred  in  the  transonic  Mach  number  region  from  Ma  =  0.80  to  0.95; 
the  largest  effects  were  evident  at  Mach  numbers  of  0.85  and  0.90.  Figure  8  presents  the 
effect  of  the  direction  of  model  sideslip  (increasing  or  decreasing  angle)  on  the  wing  upper 
surface  pressure  distributions  of  the  clean  F-4C  at  a  =  15  deg  and  0  =  2  deg  for  various 
Mach  numbers.  For  increasing  0,  the  initial  angle  was  -12  deg;  for  decreasing  /3,  the  initial 
angle  was  12  deg.  At  Mach  numbers  of  0.85  or  less,  the  flow  on  the  outboard  section  of  the 
left  wing  was  attached  on  the  forward  half  of  the  wing.  At  Ma  =  0.85,  the  rolling-moment 
coefficient  was  dependent  on  the  direction  of  model  movement  (see  Fig.  7c),  while  the 
outboard  wing  pressure  distributions  for  increasing  and  decreasing  sideslip  angle  were  Very 
similar  (see  Fig.  8c),  indicating  that  the  flow  on  the  right  wing  upper  surface  was  separated 
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with  decreasing  sideslip  angle.  At  Mach  numbers  of  0.9  and  0.95,  the  flow  was  separated  on 
the  outboard  section  of  the  left  wing  with  increasing  sideslip  angle.  With  decreasing  sideslip 
angle  (decreasing  from  0  =  12  deg),  the  flow  on  the  outboard  section  of  the  left  wing  was 
partially  reattached  on  the  forward  half  of  the  wing  (see  Fig.  8d).  Although  the  flow  was 
separated,  the  wing  pressure  distributions  at  =  0.9  and  0.95  do  not  indicate  the  presence 
of  a  shock  system  on  the  wing  upper  surface. 

4.3  DATA  REPEATABILITY 

Repeatability  of  the  aerodynamic  hysteresis  on  the  F-4C  model,  including  force  data  and 
wing  pressure  distribution  data,  indicates  that  the  phenomenon  is  not  random.  Figure  9  gives 
an  example  of  the  wing  pressure  distribution  repeatability  at  =  0.9,  a  =  15  deg,  and 
j8  =  0  deg  with  increasing  sideslip  angle  from  the  maximum  negative  value  03  =  -12  deg). 

The  effect  of  initial  model  sideslip  direction  on  the  hysteresis  characteristics  for  the  clean 
F-4C,  at  Ma  =  0.9,  a  =  15  deg,  is  shown  in  Fig.  10.  The  model  was  positioned  at  0  deg 
angle  of  attack  and  angle  of  sideslip  at  the  start  of  each  polar.  The  model  was  then  pitched  to 
a  =  15  deg  at  the  start  of  each  yaw  polar,  and  sideslip  angle  was  either  increased  or 
decreased  from  0  =  0.  Figure  10a  presents  the  effect  of  the  direction  of  model  movement  on 
the  rolling-moment  coefficient.  The  rolling-moment  coefficient  was  similar,  regardless  of 
initial  model  direction.  Figure  10b  presents  wing  upper  surface  pressure  distributions  at  0  = 
2  deg  with  a  positive  initial  model  direction.  The  wing  pressure  distributions  for  increasing 
(from  0  =  -12  deg)  or  decreasing  (from  0  =  12  deg)  0  were  the  same,  regardless  of  initial 
model  direction  when  hysteresis  was  present  (see  Fig.  8d). 

4.4  TIME  DEPENDENCE 

The  time  dependence  of  the  aerodynamic  hysteresis  phenomenon  of  the  F-4C  model  was 
also  investigated.  Figure  1 1  presents  wing  pressure  distribution  data  taken  at  selected  time 
intervals  with  tunnel  conditions  and  model  attitude  held  constant.  These  conditions  were 
Mw  -  0.9,  a  -  15  deg,  and  0  =  0  deg,  with  0  increasing  from  the  maximum  negative  value 
03  =  -12  deg).  The  model  was  initially  positioned  at  a  =  15  deg,  0  =  -12  deg.  Sideslip  angle 
was  then  increased  in  increments  of  two  degrees  as  in  a  normal  yaw  polar.  When  the  value  of 
sideslip  angle  became  zero,  the  model  attitude  was  held  constant  and  data  were  acquired 
over  a  10-min  time  interval.  The  measured  rolling-moment  coefficient  was  constant  at  Cp  « 
-0.005,  The  wing  pressure  distribution  did  not  change  significantly  with  time,  and  the  flow  at 
the  outboard  station  remained  separated.  For  the  above  test  condition  and  model  attitude 
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history,  the  hysteresis  phenomena  was  not  time  dependent.  However,  the  flow  may  be  time 
dependent  at  other  angles  of  sideslip  where  the  flow  is  either  separating  or  reattaching  (e.g., 
Fig.  10a,  j3  =  ±6  deg). 

4.5  FLOW  VISUALIZATION 

In  Fig.  12,  flow  visualization  photographs  of  the  clean  F-4C  model  show  the  effect  of  the 
previous  model  attitude  history  on  the  wing  flow  patterns  during  yaw  polars  at  Mw  =  0.9. 
For  the  oil-flow  and  tuft  tests,  the  model  was  positioned  at  a  =  15  deg  and  (3  =  -12  deg  (or 
/?  =  12  deg)  at  the  beginning  of  the  yaw  polar.  Sideslip  angle  was  increased  (or  decreased)  to 
(3  =  12  deg  (or  (3  =  -12  deg)  and  then  returned  to  =  -12  deg  (or  0  =  12  deg).  Photographs 
were  taken  at  selected  angles.  Figures  12a  and  b  present  oil-flow  and  tuft  photographs,  taken 
at  (3  =  2  deg,  for  increasing  and  decreasing  sideslip  angle,  respectively.  The  oil-flow  and  tuft 
photographs  in  Fig.  12a  indicate  flow  separation  on  the  outboard  left  wing.  This  is 
confirmed  by  the  left  wing  pressure  distribution  presented  in  Fig.  8d.  The  flow  patterns  of 
Fig.  12b  and  the  wing  pressure  distributions  of  Fig.  8d  show  that  with  decreasing  (3,  the  flow 
has  reattached  on  the  left  wing,  but  it  is  now  separated  on  the  outboard  right  wing. 
Examination  of  the  photographs  shows  that  asymmetric  flow  separation  initially  occurs  near 
the  leading  edge  of  the  wing  inboard  of  the  wing  leading-edge  notch.  The  hysteresis 
phenomenon  was  a  result  of  asymmetric  flow  separation,  which  was  dependent  on  the 
previous  model  attitude  history. 

4.6  EFFECT  OF  SIMULATED  LEADING-EDGE  SLATS 

Simulated  wing  leading-edge  slats  were  tested  on  the  basic  F-4  wing.  The  leading-edge 
slats  were  attached  to  the  basic  wing;  no  attempt  was  made  to  simulate  the  actual  wing 
leading  edge  with  slats  extended.  Figure  13  presents  a  comparison  of  wing  pressure 
distributions  for  the  configurations  tested  at  Mw  =  0.9  and  a  =  15  deg  with  increasing  a. 
For  all  configurations  except  the  basic  configuration,  1,  and  the  configuration  with 
outboard  slats  only,  6,  the  flow  on  the  outboard  wing  panel  was  separated  near  the  leading 
edge.  A  comparison  of  the  rolling-moment  coefficients  for  configurations  1  and  6,  at  MB  = 
0.9  and  a  =  15  deg,  is  given  in  Fig.  14.  Hysteresis  in  rolling-moment  coefficient  was  not 
present  for  Configuration  6.  Figure  15  presents  wing  upper  surface  pressure  distributions  for 
configurations  1  and  6,  with  increasing  and  decreasing  sideslip  angle,  and  verifies  that  the 
flow  remained  attached,  regardless  of  sideslip  direction,  on  the  wing  with  the  outer  wing 
panel  slats.  Adding  the  outboard  slat,  Si,  also  affected  the  wing  pressure  distribution  at  the 
inboard  station.  The  hysteresis  appears  to  be  caused  by  a  wing  leading-edge-related 
separation  that  originates  inboard  of  the  wing  leading-edge  notch  and  that  affects  the  flow 
on  the  wing  outer  panel. 
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5.0  MOTION  SIMULATION  STUDY  DESCRIPTION 


5.1  GENERAL 

The  importance  of  nonlinearities  in  predicting  aircraft  motion  in  the  high  angle-of-attack 
regime  has  long  been  recognized.  Reference  5  presents  the  results  of  a  study  on  motion 
prediction  of  the  F-5E  aircraft  including  the  effects  of  nonlinearities  and  aerodynamic 
hysteresis  phenomena.  Recently,  symposiums  on  high  angle-of-attack  aerodynamics  and  its 
effect  on  aircraft  motion  prediction  have  been  held  (Refs.  6,  7,  and  8).  The  effect  of 
aerodynamic  hysteresis  on  aircraft  motion  prediction,  including  wing  rock,  has  also  been 
investigated  (Ref.  9). 

The  motion  simulation  study  described  herein  was  conducted  to  assess  the  effect  of 
hysteresis  in  the  static  rolling-moment  coefficient  on  the  predicted  motion  of  the  F-4 
aircraft.  The  rolling-moment  coefficient  was  used  because  hysteresis  was  most  predominant 
in  the  rolling-moment  coefficient.  The  data  for  all  other  static  aerodynamic  coefficients  were 
based  on  data  obtained  for  increasing  angle  of  attack  or  sideslip  angle.  The  mass,  inertia, 
and  geometric  characteristics  used  in  this  investigation  for  the  full-scale  aircraft  are  listed  in 
Table  2. 

5.2  AERODYNAMIC  DATA 

The  present  study  utilized  measured  static  force  and  moment  aerodynamic  data  obtained 
from  Refs.  2,  3,  and  10.  Control  effectiveness  data,  static  stability  derivatives,  and  damping 
derivatives  were  obtained  from  Ref.  1 1 ;  all  aerodynamic  data  were  based  on  a  rigid  airplane. 
The  effects  of  angle  of  attack  and  Mach  number  were  also  included  in  these  data. 
Aerodynamic  data  were  input  in  table  look-up  form  and  linear  interpolation  was  used. 
Incremental  values  of  rolling-moment  coefficient  attributable  to  aerodynamic  hysteresis, 
dkC^|,  as  a  function  of  angle  of  attack  and  Mach  number  based  on  wind  tunnel  data,  were 
also  tabulated  (see  Table  3).  Model  movement  studies  were  conducted  to  identify  the  sideslip 
angles  at  which  hysteresis  was  present.  Separate  yaw  polars  were  run  in  either  the  positive  or 
negative  direction  of  sideslip  angle  at  a  constant  angle  of  attack.  In  either  case,  the  model 
attitude  history  was  consistent.  At  the  start  of  each  polar,  the  model  was  positioned  at  0  deg 
angles  of  attack  and  sideslip.  The  model  was  then  pitched  to  a  =  15  deg  and  sideslip  angle 
was  either  increased  or  decreased  from  j3  =  0.  Figure  16  illustrates  how  the  rolling-moment 
coefficient  was  affected  by  the  magnitude  of  sideslip  angle  when  the  model  was  returned  to 
zero  sideslip  at  Mw  =  0.9  and  a  =  15  deg.  If  the  sideslip  angle  exceeded  approximately  ±  5 
deg,  then  hysteresis  was  present,  and  the  rolling-moment  coefficient  was  shown  to  be 
dependent  on  the  direction  of  model  movement.  Also,  the  model  was  initially  positioned  at 
a  =  15  deg  and  0=12  deg  with  =  0,9;  sideslip  angle  was  decreased  to  -4  deg  and  then 
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increased  to  4  deg.  The  initial  position  of  the  model  caused  the  rolling-moment  coefficient  to 
be  on  a  hysteresis  loop.  The  rolling-moment  coefficient  remained  on  the  hysteresis  loop 
when  sideslip  angle  was  increased  to  0  =  4  deg  from  0  =  -4  deg,  as  shown  in  Fig.  17. 

5.3  METHOD  OF  ANALYSIS 

The  differential  equations  of  motion  are  presented  in  Section  1 .0  of  the  appendix.  Also 
presented  are  equations  defining  parameters  required  by  the  six-degree-of-freedom 
equations.  Two  logics  (LOGIC  1  and  LOGIC  2)  were  formulated  to  characterize  the  rolling- 
moment  coefficient  data  based  on  results  of  model  movement  studies  and  published  motion 
simulation  studies  (Refs.  5  and  9).  The  model  movement  studies  showed  that  the  hysteresis 
characteristics  of  the  static  rolling-moment  coefficient  were  such  that  the  resulting  AC^  was 
in  a  counterclockwise  sense. 

LOGIC  1  was  used  to  determine  whether  hysteresis  was  present,  based  on  the  magnitude 
of  the  aircraft  sideslip  angle.  Figure  18a,  with  a  counterclockwise  path,  illustrates  how  AQH 
was  applied  for  LOGIC  1  with  increasing  0  from  the  initial  trimmed  condition  at  point  1.  If 
0  exceeded  the  positive  value  of  0m  (point  2),  AC^  was  added  to  the  basic  rolling-moment 
coefficient  (point  3)  and  remained  in  effect  until  the  value  of  0  decreased  to  less  than  the 
negative  value  of  0m  (point  4);  AC^  was  subtracted  from  the  basic  rolling-moment 
coefficient  (point  5)  and  remained  in  effect  until  the  value  of  0  increased  to  a  value  greater 
than  the  positive  value  of  0m  (point  6).  The  resulting  path  is  hereafter  referred  to  as  LOGIC 
1,  counterclockwise  path.  Similarly,  for  decreasing  0  (LOGIC  1,  counterclockwise  path),  the 
resulting  path  would  be  the  point  sequence  1-2- 5-6- 3-4. 

LOGIC  2  was  used  to  determine  whether  an  increment  of  rolling-moment  coefficient 
attributable  to  hysteresis  was  present,  based  on  whether  the  magnitude  of  0m  was  exceeded 
and  whether  its  sign  was  dependent  on  the  sign  of  the  0  term.  Figure  18b  shows  how  AC^ 
was  applied  to  the  basic  rolling-moment  coefficient,  thus  resulting  in  a  counterclockwise 
path.  For  LOGIC  2,  if  the  magnitude  of  0m  was  exceeded,  then  ACfH  was  applied  through 
the  maneuver,  and  the  value  of  the  rolling-moment  coefficient  jumped  directly  from  the 
increasing  0  path  to  the  decreasing  0  path  as  0  changed  sign. 

Additional  logics  were  developed  to  apply  AC^  clockwise  for  LOGIC  1  and  LOGIC  2  as 
shown  in  Figs.  18a  and  b,  respectively. 

For  the  analytical  simulations,  the  initial  flight  conditions  were  specified  and  the  control 
deflections  and  thrust  were  determined  to  obtain  trimmed  flight.  A  control  surface  was  then 
moved  to  begin  the  maneuver.  The  computer  program  was  numerically  integrated  with  time 
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by  using  a  Runge-Kutta  integration  algorithm  (Ref.  12)  to  predict  a  few  initial  points.  After 
these  initial  points  were  generated,  an  Adams-Moulton  predictor-corrector  integration 
algorithm  was  used  for  the  numerical  integration  at  each  subsequent  time  step.  Both 
algorithms  used  a  fixed  time  step  of  0.01  sec.  The  computer  program  was  written  so  that  any 
of  the  logics  described  previously  could  be  applied  with  the  option  of  the  absence  of 
hysteresis  effects.  An  examination  of  the  computed  motions  for  each  logic,  when  compared 
to  the  case  without  hysteresis  effects,  indicated  the  effects  of  hysteresis. 

6,0  STUDY  RESULTS  AND  DISCUSSION 

The  primary  area  of  investigation  was  aircraft  maneuvers  at  transonic  Mach  numbers. 
Because  aerodynamic  hysteresis  was  most  predominant  at  =  0.9,  this  Mach  number  was 
used  as  one  of  the  initial  trimmed  flight  conditions  for  the  maneuvers  described  in  the 
following  sections.  Other  flight  conditions  included  level  and  3-g  turning  flight.  Calculated 
time  histories  were  plotted  with  use  of  a  0.25-sec  time  interval. 

6.1  GENERAL 

Flight  motion  simulation  studies  were  performed  with  use  of  the  F-4C  aircraft 
characteristics  given  in  Table  2.  The  thrust  and  control  surface  deflections  required  for 
trimmed  flight  at  an  altitude  of  30,000  ft  and  at  a  Mach  number  of  0.9  were  obtained  by 
using  the  aerodynamic  data  described  in  Section  5.2.  For  the  3-g  turning  flight  condition,  it 
was  assumed  that  trimmed  thrust  could  be  obtained. 

Control  surface  deflections  for  the  pull-up  maneuver  were  a  stabilator  ramp  of  -3  deg 
from  the  trimmed  stabilator  setting  at  -1  deg/sec  from  time  equal  0  to  3  sec.  The  stabilator 
was  then  held  constant  for  the  duration  of  the  maneuver.  A  rudder  doublet  of  ±  10  deg  with 
a  rudder  deflection  rate  of  ±  20  deg/sec  was  then  performed  to  excite  the  motion  in  the 
lateral  plane.  For  turning  flight,  the  control  surface  perturbation  was  a  rudder  doublet  of 
±  10  deg  from  the  trim  rudder  setting  with  a  deflection  rate  of  ±  20  deg/sec  to  excite  motion 
in  the  lateral  plane. 

The  resulting  time  histories  of  the  aircraft  motion  were  analyzed  to  assess  the  effect  of 
hysteresis  in  the  static  rolling-moment  coefficient  on  predicted  motion.  A  value  of  |  /3m  I  =5 
was  used  for  all  analyses  except  the  j3m  study. 

6.2  PULL-UP  MANEUVER 

Figure  19  compares  the  effects  of  hysteresis  logic  on  the  aircraft  motion  (including  the 
case  without  hysteresis)  during  a  pull-up  maneuver  with  a  rudder  doublet.  Figure  19a  shows 
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a  time  history  of  Mach  number,  and  Figs.  19b,  c,  and  d  present  angle  of  attack,  sideslip 
angle,  and  roll  angle,  respectively.  As  expected,  the  aircraft  slows  down  and  angle  of  attack 
is  increased  for  all  logics.  The  application  of  LOGIC  1,  with  a  counterclockwise  AC,^,  does 
not  result  in  any  significant  change  in  the  sideslip  angle  time  history,  compared  to  the  case 
without  hysteresis,  but  there  is  significant  change  in  the  roll  angle  time  history  (see  Fig.  19d). 
This  change  in  roll  angle  time  history  results  from  the  hysteresis  which  is  present  from 
time  =5  to  ~  10  sec,  when  Mach  number  becomes  less  than  0.7  (/JM  <  -5  deg,  a  >  5  deg, 
M#  >  0.7;  see  Table  5).  The  roll-off  in  the  negative  direction  is  a  result  of  being  on  the 
lower  side  of  the  hysteresis  loop. 

The  use  of  LOGIC  1  with  a  clockwise  direction  results  in  an  increase  of  sideslip  and  roll 
angle,  compared  to  the  case  without  hysteresis,  between  time  equals  6  and  12  sec;  this 
increase  results  from  (3’s  exceeding  ±5  deg,  which  results  in  switching  sides  of  the  hysteresis 
loop.  The  roll-off  in  the  positive  direction  results  from  being  on  the  upper  side  of  the 
hysteresis  loop  until  Mach  number  decreases  to  less  than  0.7,  at  which  point  no  hysteresis  is 
present. 

Applying  LOGIC  2  with  a  counterclockwise  direction  results  in  a  decrease  in  sideslip 
angle  oscillation  and  a  change  in  the  roll  angle  time  history,  compared  to  the  case  without 
hysteresis.  This  clockwise  application  of  LOGIC  2  results  in  an  increase  in  sideslip  and  roll 
angles,  compared  to  the  case  having  no  hysteresis.  In  addition,  the  roll  angle  oscillation 
between  lime  equal  6  and  1 1  sec  increases  because  AQH  changes  signs  as  0  changes  sign.  The 
sideslip  angle  oscillates  until  Mach  number  decreases  to  below  0.7,  at  which  hysteresis  is 
absent  in  the  rolling-moment  coefficient. 

6.3  TURNING  FLIGHT 

Flight  motions  were  calculated  by  applying  a  rudder  doublet  for  the  aircraft  initially 
trimmed  in  a  3-g  turn.  Figure  20  gives  a  comparison  of  the  effects  of  how  AC^  is  applied 
with  LOGIC  1  on  the  predicted  aircraft  motion.  Traces  of  the  calculated  Mach  number,  the 
angle  of  attack,  and  the  sideslip  and  roll  angles  are  presented  for  both  the  cases  of 
counterclockwise,  and  clockwise  applications  of  AC^  and  the  case  without  hysteresis.  The 
Mach  number  and  angle  of  attack  for  the  first  6  sec  are  similar.  After  6  sec,  for  the 
counterclockwise  case,  the  Mach  number  decreased  more  rapidly  in  comparison  to  the  case 
without  hysteresis,  while  angle  of  attack  increased  until  time  equaled  approximately  15  sec. 
Beyond  15  sec,  the  angle  of  attack  decreased.  With  clockwise  application  of  AC^,  the  Mach 
number  again  decreased  more  rapidly  and  the  angle  of  attack  continued  to  increase  to  an 
angle  for  which  accurate  data  were  not  available,  7,  and  the  calculation  was  discontinued. 
The  overall  sideslip  and  roll  angle  characteristics  are  changed  significantly.  The  roll-off  in 
the  negative  direction  is  a  result  of  being  on  the  lower  side  of  the  hysteresis  loop. 
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The  wind  tunnel  test  results  indicate  that  the  hysteresis  characteristics  of  the  rolling- 
moment  coefficient  were  such  that  the  coefficient  tended  to  remain  on  the  upper  or  lower 
portion  of  the  loop,  depending  on  the  direction  of  model  movement  and  the  magnitude  of 
the  sideslip  angle  (see  Figs,  16  and  17).  Based  on  the  wind  tunnel  results  (discussed  in  Section 
4.7)  a  value  of  |  0m  I  =  5  deg  was  used  for  the  basic  study.  The  resulting  calculated  motion 
with  LOGIC  1  is  a  roll-off,  with  the  direction  dependent  upon  on  which  side  of  the  hysteresis 
loop  the  simulation  ended.  Flight  motions  were  calculated  for  a  reduced  value  of  0m  with 
LOGIC  1.  Figure  21  compares  the  calculated  motions  with  the  counterclockwise  and 
clockwise  applications  of  AC^  for  |  0M  |  =  5  and  2  deg.  For  the  counterclockwise  case  (Fig. 
21a),  there  are  significant  differences  in  all  of  the  parameters  (M0,  a,  0,  <p),  with  the  roll 
angle  dependent  upon  on  which  side  of  the  hysteresis  loop  the  simulation  ended.  In  the 
clockwise  case  (Fig.  2  lb)  with  |  0m  |  =  2  deg,  the  presence  of  hysteresis  caused  oscillations  in 
sideslip  and  roll  angle,  while  for  the  clockwise  application  with  |  0M  |  =5  deg,  the  hysteresis 
caused  a  roll-off  in  the  negative  direction.  A  comparison  of  the  calculated  motions  using 
LOGIC  I  (with  both  counterclockwise  and  clockwise  application  of  AC^)  and  the  effect  of 
the  magnitude  of  0m  clearly  demonstrates  the  importance  of  how  static  aerodynamic 
hysteresis  is  interpreted  and  applied  in  motion  simulation  calculations.  It  should  be 
remembered  that  the  static  wind  tunnel  data  show  that  hysteresis  was  present  with  a 
counterclockwise  direction  and  that  the  magnitude  of  |  0m  j  was  approximately  5  deg  (see 
Fig.  10). 

LOGIC  2  is  used  in  additional  flight  motion  calculations  in  which  the  sign  of  AQy,  is 
dependent  on  the  sign  of  0  once  |  0M  |  has  been  exceeded.  A  comparison  of  the  effect  of 
how  ACfn  is  applied  with  LOGIC  2  on  the  predicted  aircraft  motion  is  given  in  Fig.  22.  A 
comparison  of  the  case  without  hysteresis  and  the  case  of  counterclockwise  application  of 
AQii  reveals  the  following:  For  the  case  with  hysteresis  present,  the  Mach  number  decreased 
more  rapidly  as  the  angle  of  attack  tended  to  increase  more  rapidly  to  approximately  12  deg 
at  time  equal  to  11  seconds.  There  were  significant  differences  in  the  sideslip  and  roll 
damping,  as  hysteresis  caused  increased  damping.  Comparing  the  case  without  hysteresis 
and  the  case  of  the  clockwise  application  of  AC^  shows  the  following:  For  the  case  with 
hysteresis,  the  Mach  number  tended  to  remain  near  the  initial  condition  (M#  =  0.9),  while 
the  angle  of  attack,  sideslip  angle,  and  roll  angle  oscillated  with  significantly  reduced 
damping. 

6.4  MOTION  SENSITIVITY 

The  analysis  described  in  the  previous  sections  was  based  on  a  six-degree-of-freedom 
motion  simulation  computer  program  that  was  modified  to  include  aerodynamic  hysteresis 
in  the  static  rolling-moment  coefficient.  The  initial  analysis  was  made  with  use  of  LOGIC  1 
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and  LOGIC  2  with  a  counterclockwise  application  of  AC^.  This  counterclockwise  direction 
was  based  on  static  wind  tunnel  test  data;  the  application  of  AQH  resulted  in  increased  yaw 
and  roll  damping.  Clockwise  application  of  ACtH  resulted  in  hysteresis-induced  motion. 

7.0  CONCLUDING  REMARKS 

In  this  investigation  of  the  effects  of  the  hysteresis  phenomenon  with  the  F-4,  the 
following  were  observed: 

1.  The  results  of  wind  tunnel  tests  with  a  0.05-scale  F-4  model  showed  that 
aerodynamic  hysteresis  was  present  in  the  static  longitudinal  and 
lateral/directional  data.  Hysteresis  occurred  primarily  in  the  rolling-moment 
coefficient  but  was  also  present  in  the  pitching-moment  and  normal-force 
coefficients.  The  static  aerodynamic  characteristics  were  dependent  on  the 
model  attitude  history  and  were  repeatable  for  identical  a/0  sweeps  in  the  wind 
tunnel. 

2.  The  hysteresis  phenomenon  exhibited  in  the  aerodynamic  data  was  reflected  in 
the  flow  over  the  wing  of  the  F-4  model. 

3.  Flight  motion  simulations  in  which  the  hysteresis  characteristics  of  the  rolling- 
moment  coefficient  were  modeled  on  the  direction  of  the  static  wind  tunnel  data 
produced  indications  of  increases  in  yaw  and  roll  damping. 

4.  Flight  motion  simulations  with  the  hysteresis  characteristics  of  the  rolling- 
moment  coefficient  modeled  against  the  direction  of  the  static  wind  tunnel  data 
resulted  in  hysteresis-induced  oscillating  motion  with  decreases  in  yaw  and  roll 
damping. 

5.  Based  on  the  preceding  two  observations,  it  can  be  concluded  that  flight 
simulations  using  limited  static  wind  tunnel  data  can  be  misleading,  particularly 
in  regard  to  hysteresis  effects  associated  with  the  rolling-moment  coefficient. 
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a.  General  arrangement 

Figure  1.  Details  and  dimensions  of  the  F-4C  model. 
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Figure  1.  Concluded. 
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a.  F-4C,  Configuration  1 
Figure  2.  Model  installation  in  Tunnel  4T. 
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c.  F-4C,  Configuration  6 
Figure  2.  Concluded. 
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b.  Pitching-moment  coefficient 
Figure  3.  Continued. 
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d.  Side-force  coefficient 
Figure  3.  Continued. 


30 


AEDC-TR-80-10 


Sym 

O  0.70 

n  0.80 

<£k  0.85 

O  0.90 
y  0.95 


f.  Rolling-moment  coefficient 
Figure  3.  Concluded. 
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b.  a  =  8  to  15  deg 
Figure  4.  Continued. 
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c.  a  =  IS  to  20  deg 
Figure  4.  Concluded. 
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b.  Mm  =  0.8,  a  =  15  deg 
Figure  5.  Continued. 
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d.  =  0.9,  a  =  15  deg 
Figure  5.  Continued. 
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Figure  5.  Concluded. 
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Figure  8.  Effect  of  model  sideslip  direction  on  wing  upper 

surface  pressure  distributions,  a  =  15  deg,  0  =  2  deg. 
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Figure  8.  Continued. 
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c.  =  0.85 
Figure  8.  Continued. 
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Figure  8.  Concluded. 
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Figure  9.  Wing  upper  surface  pressure  distribution  repeatability, 
M„  =  0.9,  a  =  15  deg,  (3  =  0  deg. 
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b.  Increasing  f?,  from  maximum  negative  angle 
Figure  9.  Concluded. 
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hysteresis  characteristics,  M„  =  0.9,  o  =  15  deg. 
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Fiaure  10.  Concluded. 
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Figure  11.  Wing  upper  surface  pressure  distribution  time 
dependence,  with  increasing  sideslip  angle, 

=  0.9,  a  -  15  deg,  (3  =  0  deg. 
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a.  Increasing  |3  b.  Decreasing  (i 

Figure  12.  Effect  of  previous  model  attitude  history 

on  the  wing  flow  patterns  of  the  clean  F-4C 
=  0.9,  a  =  15  deg,  (3=2  deg. 
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Figure  13.  Effect  of  slats  on  wing  upper  surface 
pressure  distributions  with  increasing 
angle  of  attack,  =  0.9,  a  =  15  deg. 
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□  6  (F-4C  with  Outboard  Slat,  S^) 


6 

Figure  14.  Effect  of  outboard  slat,  Si ,  on  the  rolling-moment 
coefficient,  =  0.9,  a  -  15  deg. 
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Figure  15.  Effect  of  outboard  slat,  Si ,  on  wing  upper  surface 
pressure  distribution,  =  0.9,  a  =  15  deg, 

/?  =  2  deg. 
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lure  16.  Effect  of  the  magnitude  of  sideslip  angle  and  the 

direction  of  model  movement  on  the  rolling-moment 
coefficient.  M  =  0.9,  a  =  15  deg. 
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Figure  17.  Characteristics  of  the  rolling-moment  coefficient 
in  a  hysteresis  loop,  M,,  =  0.9,  a  =  15  deg. 
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Counterclockwise  Path 


b.  LOGIC  2 

Figure  18.  Motion  simulation  hysteresis  logics. 
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Time 


a.  Mach  number 

Figure  19.  The  effect  of  hysteresis  logic  on  the  aircraft 
motion— pull-up  maneuver  from  level  flight. 


0  2  4  6  8  10  12  14 

Time 


b.  Concluded 
Figure  19.  Continued. 
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c.  Sideslip  angle 
Figure  19.  Continued. 
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c.  Concluded 
Figure  19.  Continued. 
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d.  Roll  angle 
Figure  19.  Continued. 
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d.  Concluded 
Figure  19.  Concluded. 
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b.  Angle  of  attack 
Figure  20.  Continued. 
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d.  Roll  angle 
Figure  20.  Concluded. 
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b.  Clockwise  ACg^ 
Figure  21.  Concluded. 
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Figure  22.  Effect  of  ACg^  direction  on  the  aircraft 
motion  with  LOGIC  2,  3-g  turning  flight. 
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Figure  22.  Continued. 
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Time 


c.  Sideslip  angle 
Figure  22.  Continued. 
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d.  Roll  angle 
Figure  22.  Concluded. 
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Table  1.  Configuration  Identification 

Conf igurat ion 

Number  Configuration 

1  Clean  F-4C  model 

2  F-4C  with  S1  and 

3  F-4C  with  S3 

4  F-4C  with  (not  tested  -  see  Ref.  2) 

5  F-4C  with  and  S£ 

6  F-4C  with 
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Table  2.  Mass,  Inertia,  and  Dimensional  Characteristics 
of  the  Full-Scale  F-4  Aircraft 


Load ing 


Weight,  lb.  36,950 

Cg,  percent  Mean  Aerodynamic  Chord  33.0 


2 

Moments  of  Inertia,  Slug-ft 

Ix  22,600 

IY  163,100 

Iz  182,000 

:xz  5’450 

Wing 

Span,  ft 
Area,  ft2 

Mean  Aerodynamic  Chord,  ft 


38.41 

538.34 

16.04 


Inertia  Coupling  Terms 


XY~ 


-0.836 


0.241 


0.977 


-0.772 


0.033 

0.030 
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Table  3.  Incremental  Rolling-Moment  Coefficient 
Attributable  to  Hysteresis 


Mach 

Angle  of  attack,  deg 

Number 

0 

5 

10 

15 

20 

0.7 

0 

0 

0 

0 

0 

0.8 

0 

0 

0 

0.005 

0 

0.85 

0 

0 

0.001 

0.009 

0 

0.9 

0 

0 

0.002 

0.006 

0 

0.95 

0 

0 

0.001 

0.005 

0 

1.1 

0 

0 

0 

0.002 

0 
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APPENDIX  A 
EQUATIONS  OF  MOTION 

A-1.0  S1X-DEGREE-OF-FREEDOM, 

RIGID-BODY  EQUATION  OF  MOTIONS 

The  dynamic  equations  required  to  specify  the  translational  and  rotational  motions  of  a 
rigid  body  moving  through  space  are  described  in  this  appendix.  The  six-degree-of-freedom 
nonlinear  differential  equations  representing  the  linear  and  angular  accelerations  on  a 
moving  body  axis  system  (Fig.  A-2)  whose  origin  is  at  the  aircraft  center  of  mass  are  given 
below. 

Forces: 

f  tx 

u  =  rv  —  qw  -  a  sin  6  +  - —  T  — 

m  m 

F 

v  =  pw  -  ru  +  g  cos#  sin<£  + 

F  T 

m  /i  j  *  * 

w  =  qu  —  pv  —  e  cos  tr  cos  0  +  —  +  — 
r  r  m  m 

Moments: 


’y  ~  *2 


*xz  x  Mx 

P  =  — - -  qr  -*■  — —  (r  +  pq)  — 

X  x  X 


.  lz  !X  *XZ  ,  ,  My  MYT 

q  =  ~i7~pr  +  v  ~p  )  +  V  — 


‘x  _  *Y  ‘xz  1V,Z 

r  =  — -  pq  +  -jj —  (  P  -  qr)  +  — 

!Z  'z  2 

The  external  forces  and  moments  (Fx,  FY,  Fz,  Mx,  MY,  and  Mz)  in  the  equations  are 
comprised  of  aerodynamic  coefficients  for  the  aircraft.  The  external  force  and  moment 
contributions  attributable  to  engine  thrust  are  represented  by  Tx,  TZ(  and  MYT  and  are 
developed  in  Section  A-2.0.  Development  of  the  aerodynamic  math  model  used  is  presented 
in  Section  A-3.0. 
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The  auxiliary  equations  used  in  this  analysis  are  the  following: 


a  = 


p  -  (f) 


a 


u  w  —  wu 

2  2 
u  +  w 


V 


(Ull  +  vv  +  ww  \ 


V  =  yju2  +  v2  +  w2  , 
6  =  q  cos  <fi  —  t  sin  <f> 


y  = 


p  +  tan  6  (r  cos  <f>  +  q  sin  <f>) 

T  cos  0  +  q  ain  <f> 
cos  6 


The  equations  of  motion  are  numerically  integrated  to  provide  time  histories  of  the  aircraft 
motion. 


A-2.0  EQUATIONS  DEFINING  THE  FORCE  AND  MOMENT 
CONTRIBUTIONS  ATTRIBUTABLE  TO  ENGINE  THRUST 


deg 

Line 


Figure  A-1.  External  force  and  moment  contributions 
attributable  to  engine  thrust. 
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Forces 


Moments 


Tx  =  T  cos  e 


Mx-p  —  0 


Ty  =  0  Myt  =  T(d) 

T2  =  -T  sin  ( 


where 

d  =  Distance  between  thrust  line  and  CG 

Assumptions 

1.  Engine  thrust  line  is  parallel  to  X-Z  plane. 

2.  There  is  no  variation  in  Y  and  Z  center  of  gravity  location. 

A-3.0  EQUATIONS  DEFINING  THE  TOTAL  AERODYNAMIC 
DATA  ALONG  AND  ABOUT  EACH  BODY  AXIS 

Longitudinal  Axis  Plane 

FX  =  ‘loa5  I  C\r(a>  0*  +  ACx{a,  Mm,  8h)  +  f"cx  (a)"j  — 

i  L  n  J  2V 

FX  =  <US  Cz(0,/?,  MJ  +  ACz(a,MaotSH)+  fcz  (MJq  +  Cz.(MJil^. 

1  U  q  J  2V 

MY  =  <USC  |  C>.  &  MJ  +  ACm(a,  Mm,  5h)  +  [Cm?(Mj  q  +  Cm.  <MJ  d]  ~  | 

Lateral-Directional  Axis  Plane 


+ 


Pb 
2  v 


+ 


[ 


C  (a,  M  ) 

00 ' 
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Mz  =  q-Sb  |  cn(«,  0,  MJ  +  AC>,  Mw,  fir)  +  ACn(a,  Mm.  SJ 

*[c»,la,MJ]  « +  [c»>"-)]iv| 

Mx  =  qJSb  |  Cp(a,  /3,  Mn)  +  AC^(a,  JS,  SJ  +  ACg(a,  Sr) 

+  ACjMa,  MJ  +  [c^(a.  MjJ  ^  -  [^Gg  (a,  Ms)j  ~  j 


Center  of  Gravity 


Z  1  V 

Figure  A-2.  Body  axial  system. 
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NOMENCLATURE 

BL  Buttline,  in. 

b  Wing  span;  model  =  1.927  ft,  full  scale  =  38.41  ft 

Ca  Axial-force  coefficient,  measured  axial-force/q^S 

CG  Center-of-gravity  location,  percent  chord 

Q.  Centerline 

Q  Rolling-moment  coefficient,  rolling  moment/qoaSb 

Qp  Derivative  of  rolling-moment  coefficient  with  respect  to  roll  rate, 

3Cf/d(pb/2V),  per  radian 

Cfr  Derivative  of  rolling-moment  coefficient  with  respect  to  yaw  rate, 

3Q/d(rb/2V),  per  radian 

Cm  Pitching-moment  coefficient,  pitching  moment/qwSb 

Cmq  Derivative  of  pitching-moment  coefficient  with  respect  to  pitch  rate, 

flCm/3(qc/2V),  per  radian 

Cm(j  Derivative  of  pitching  moment  with  respect  to  a,  3Cm/3(dc/2V),  per  radian 

Cn.Cz  Normal-force  coefficient,  (CN  =  -Cz)  normal  force/q^S 

Cn  Yawing-moment  coefficient,  yawing  moment/q^Sb 

C„p  Derivative  of  yawing-moment  coefficient  with  respect  to  roll  rate, 

3Cn/3(pb/2V),  per  radian 

C„r  Derivative  of  yawing-moment  coefficient  with  respect  to  yaw  rate, 

3Cn/3(rb/2V),  per  radian 

Cp  Pressure  coefficient,  p  -  P^/q^ 

Cx  Longitudinal-force  coefficient,  longitudinal  force/q^S 

Cxq  Derivative  of  longitudinal-force  coefficient  with  respect  to  pitch  rate, 

3Cx/3(qc/2V),  per  radian 

Cy  Side-force  coefficient,  side  force/q^S 

CyP  Derivative  of  side-force  coefficient  with  respect  to  roll  rate,  dCY/d(pb/2V), 

per  radian 
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Derivative  of  side-force  coefficient  with  respect  to  yaw  rate,  dCy/d(rb/2V), 
per  radian 

Derivative  of  normal-force  coefficient  with  respect  to  pitch  rate, 
9Cz/d(qc/2V),  per  radian 

Derivative  of  normal-force  coefficient  with  respect  to  a,  3Cz/9(ac/2V), 
per  radian 

Local  wing  chord 

Mean  aerodynamic  chord,  scale  model  =  0.802  ft,  full  scale  =  16.04  ft 
Distance  from  thrust  line  to  center  of  gravity,  ft 
Fuselage  station,  in. 

Force  acting  along  X-body  axis,  lb 
Force  acting  along  Y-body  axis,  lb 
Force  acting  along  Z-body  axis,  lb 
Acceleration  due  to  gravity,  ft/sec2 
Altitude,  ft 

Moments  of  inertia  about  X-,  Y-,  and  Z-body  axes,  respectively;  slug-ft2 
Product  of  inertia,  slug-ft2 
Mass,  slugs 

Moment  acting  about  X-body  axis,  ft-lb 

Moment  acting  about  Y-body  axis,  ft-lb 

Moment  acting  about  Y-body  axis  caused  by  engine  thrust,  ft-lb 

Moment  acting  about  Z-body  axis,  ft-lb 

Free-stream  Mach  number 

Local  static  pressure,  psfa 

Roll,  pitch,  and  yaw  rates  about  X-,  Y-,  and  Z-body  axes,  respectively, 
radians/ sec 
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p,q,r  Roll,  pitch,  and  yaw  accelerations  about  X-,  Y-,  and  Z-body  axes, 

respectively;  radians/sec2 

pt  Free-stream  total  pressure,  psfa 

Free-stream  static  pressure,  psfa 

qB  Free-stream  dynamic  pressure,  psf 

Re  Unit  Reynolds  number,  per  foot 

S  Wing  reference  area;  scale  model  =  1.325  ft2;  full  scale  =  538.34  ft2 

T  Engine  thrust,  lb 

Tx  Component  of  thrust  along  X-body  axis,  lb 

Tz  Component  of  thrust  along  Z-body  axis,  lb 

u,v,w  Components  of  total  velocity  along  X-,  Y-,  and  Z-body  axes,  respectively; 

ft/sec 

V  Total  velocity,  ft/sec 

WL  Waterline,  in. 

X,Y,Z  Body  axes 

x,y,z  Linear  distance  along  X-,  Y-,  and  Z-body  axes,  respectively,  ft 

x/c  Fraction  of  wing  chord  measured  from  the  leading  edge 

y/c  Fraction  of  wing  chord  measured  vertically  from  chord  line 

a  Angle  of  attack  referenced  to  model  waterline,  deg 

3  Angle  of  sideslip,  deg 

3m  Angle  of  sideslip  at  which  hysteresis  is  present,  deg 

<5a  Aileron  deflection  (positive  when  trailing  edge  of  right  aileron  is  down),  deg 

6h  Stabilator  deflection  (positive  when  trailing  edge  is  down),  deg 

6r  Rudder  deflection  (positive  when  trailing  edge  is  left),  deg 

e  Angle  between  thrust  line  and  X-body  axis,  deg 
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e 

* 

SUBSCRIPT 


Angle  between  X-body  axis  and  horizontal  measured  in  vertical  plane,  deg 

Incremental  rolling-moment  coefficient  attributable  to  hysteresis, 
incremental  rolling  moment/q^Sb 

Angle  between  Y-body  axis  and  horizontal  measured  in  vertical  plane,  deg 
Angle  between  Y-body  axis  and  vertical  measured  in  a  horizontal  plane,  deg 

Derivative  with  respect  to  time 
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